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Abstract 



The inorganic crystal scintillator CsI(Tl) has been used for low energy neutrino 
and Dark Matter experiments, where the intrinsic radiopurity is an issue of major 
importance. Low-background data were taken with a CsI(Tl) crystal array at the 
Kuo-Sheng Reactor Neutrino Laboratory. The pulse shape discrimination capabil- 
ities of the crystal, as well as the temporal and spatial correlations of the events, 
provide powerful means of measuring the intrinsic radiopurity of 137 Cs as well as 
the 235 U, 238 U and 232 Th series. The event selection algorithms are described, with 
which the decay half-lives of 218 Po, 214 Po, 220 Rn, 216 Po and 212 Po were derived. 
The measurements of the contamination levels, their concentration gradients with 
the crystal growth axis, and the uniformity among different crystal samples, are 
reported. The radiopurity in the 238 U and 232 Th series are comparable to those of 
the best reported in other crystal scintillators. Significant improvements in mea- 
surement sensitivities were achieved, similar to those from dedicated massive liquid 
scintillator detector. This analysis also provides in situ measurements of the detec- 
tor performance parameters, such as spatial resolution, quenching factors, and data 
acquisition dead time. 

PACS Codes: 29.40.Mc, 29.25.R, 07.05.Kf. 

Keywords: Scintillation detector, Radioactive sources, Data analysis. 
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1 Introduction 



The merits and physics potentials for crystal scintillators in low-energy low-background 
experiments were recently discussed [1]. This detector technology has been adopted in 
double beta decay dark matter and reactor neutrino experiments, and considered for low 
energy solar neutrino detection. An important common issue is the contamination of 
radioactive isotopes in the crystals. This article reports measurements of the intrinsic 
radiopurity in CsI(Tl) crystals. 

The measurements were performed at the Kuo-Sheng (KS) Nuclear Power Plant where 
a laboratory located at a distance of 28 m from the reactor core has been built to 
study low energy neutrino physics [2]. A multi-purpose "inner target" detector space 
of 100 cmx80 cmx75 cm is enclosed by An passive shielding materials which have a to- 
tal weight of 50 tons. The shielding provides attenuation to the ambient neutron and 
gamma background, and consists of, from inside out, 5 cm of OFHC copper, 25 cm of 
boron-loaded polyethylene, 5 cm of steel, 15 cm of lead, and cosmic-ray veto scintillator 
panels. 

Different detectors can be placed in the inner space for different scientific goals. The 
detectors are read out by a versatile electronics and data acquisition (DAQ) system[3] 
based on 16-channel, 20 MHz, 8-bit Flash Analog-to-Digital-Convertor (FADC) modules. 
The readout allows full recording of all the relevant pulse shape and timing information for 
as long as 500 fxs after the initial trigger with a resolution of 1 /is. The data are then read 
out and transferred to the host computer, and the average DAQ dead time is 5.13 ms per 
event. The timing information for every event was recorded with a resolution of ~1 ms. 
The time differences between subsequent events and between different triggers within an 
event can be evaluated. Software procedures have been devised to extend the effective 
dynamic range from the nominal 8-bit measurement range provided by the FADC [4]. 

One of the detectors placed inside the inner target volume is an array of CsI(Tl) crystal 
scintillators^], the schematic of which is displayed in Figure 1. Each crystal module is 
2 kg in mass and consists of a hexagonal- shaped cross-section with 2 cm side and a length 
of 40 cm. The light output are read out at both ends (Qr and Ql) by custom-designed 
29 mm diameter photo-multiplier tubes(PMTs) with low-activity glass. The sum and 
difference of the PMT signals give information on the energy and the longitudinal Z- 
position of the events, respectively. Performance of the prototype modules was published 
in Ref. [6]. 

A total of 186 kg (or 93 modules) was commissioned for the 2003 data taking. The 
physics goal is to measure the Standard Model neutrino-electron scattering cross sections, 
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and thereby to provide a measurement of sin 2 6> w at the untested MeV range. The strat- 
egy [7] is to focus on data at high (>2 MeV) recoil energy where the uncertainties due 
to the reactor neutrino spectra are small. The large target mass compensates the drop in 
the cross-sections at high energy. 

This article focuses on the measurements of the intrinsic 137 Cs level as well as those for 
the 235 U, 238 U and 232 Th series in the CsI(Tl) crystals using the correlated events from j3-a 
and a-a decay sequences. Our work complements and extends the scope of those previous 
efforts on this [8, 9] and other crystal scintillators: Nal(Tl) [11, 12, 13], CdW0 4 [14], 
GSO(Ce) [15], CeF 3 [16], BaF 2 [17] and CaF 2 (Eu) [18]. The achieved sensitivities match 
the best among these results as well as those from massive liquid scintillator facility [19] . 
Several detector parameters can be derived from these measurements, providing a means 
of in situ monitoring of the detector calibration and performance. 

In addition, another CsI(Tl) detector array of mass 200 kg is being prepared in the 
Yang- Yang Underground Laboratory in South Korea to look for Cold Dark Matter [9]. 
Knowledge of the intrinsic radioactivity in CsI(Tl) is essential to both experiments. The 
methods devised are applicable to a wide spectrum of low-background experiments where 
a/7 discrimination is possible and the detailed timing information are available. 



2 Intrinsic 67 Cs Level 



The isotope 13 Cs is produced artificially as fission waste from power reactors and atomic 
weapon tests. Cesium compounds are highly soluble and hence 137 Cs can easily be intro- 
duced into cesium products during the processing from cesium ore (like pollucite). Once 
introduced, there is no chemical means to separate 137 Cs from the stable 133 Cs. This is 
the dominant one background contribution to Dark Matter experiments using CsI(Tl) 
crystals as target [9]. 

The isotope decays via [10]: 

137 Cs -> 137 Ba* + e~ (ri = 30.07 y) , 
137 Ba* -> 137 Ba + 7 (E 1 = 662 keV) 

with the emission of a 7-ray of energy 662 keV. A total of 31.3 kg-day of CsI(Tl) data was 
analyzed. The spectrum is displayed in Figure 2a, showing a distinct line at this energy 
with a root-mean-square (RMS) resolution of 5.1%. Other naturally-occuring 7-lines, like 
those from 40 K and 208 T1, can also be observed. 

The Z-position distribution of the 137 Cs events are shown in Figure 2b. Events within 
±3er of the 7-peak give excess on both ends of the 40-cm long crystals due to external 
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background. The uniform distribution after a ±1(7 cut indicates the source of 137 Cs is 
internal to the crystals. 

Simulation studies give an efficiency of 42% for full capture of the 662 keV 7-ray in 
a single crystal. Accordingly, an average activity of 61±2 mBq kg" 1 , or equivalently a 
contamination level of (1.55 ± 0.05) x 10" 17 g/g, can be derived for the CsI(Tl) crystals. 
The RMS spread of the activity among the 31 crystals studied is 30%. 

Alternatively, the radiopurity of the Csl powder with which the CsI(Tl) crystals were 
produced have been measured with a high-purity germanium detector. The measured 
137 Cs abundance is (1.7 ± 0.3) x 10~ 17 g/g, consistent with the in situ results. 



3 Pulse Shape Discrimination in CsI(Tl) Crystals 

It has been well-studied [6, 20] that the light emission profiles of scintillating CsI(Tl) crys- 
tals exhibit different shape for 7-rays and electrons (that is, minimum ionizing particles), 
as compared to that for a-particles and nuclear recoils at the >100 keV regime . Heavily 
ionizing events due to a-particles and nuclear recoils have faster decays than those from 
e/7's — opposite to the response in liquid scintillator [21]. The average pulse shapes for 
both categories are depicted in Figures 3, where t=0 is defined by the trigger instant. This 
characteristic property makes particle identification possible with this scintillator [22]. 

Matured pulse shape discrimination (PSD) techniques have been devised at high ener- 
gies where the photo-electrons are abundant. The "double charge method" [23] involves 
the comparison of the "total charge" (Q to t) and the "partial charge" (Qpart), which are 
the total and partial integration of the pulse, respectively. In this analysis, Q to t and Q pa rt 
were evaluated by integrating the entire pulse from —0.75 to 12.8 /is, and the tail from 5 
to 12.8 /is, respectively, with respect to the timing shown in Figure 3. 

With the digitized pulse information measured by the FADCs, the mean time method 
is often used. The mean time is defined as 

E (Ai ti) 

<*> = ^tat • (1) 

i 

where A ; is the FADC-amplitude at time-bin t;. The PSD capabilities to make a/7 
identification for MeV events are displayed in Figure 4a&b for both methods, indicating 
excellent separation. These features can be extended to the low energy regime where 
photo-electron statistics are limited, via advanced software techniques [24]. This is rele- 
vant to the Dark Matter experiments. 
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4 Intrinsic Radiopurity of U/Th Series 



The 35 U, U and 232 Th decay series are naturally occurring, and exist in all materials 
at some levels. The processing of the raw materials to CsI(Tl) crystals can preferentially 
introduce or remove certain isotopes within these series such that secular equilibrium can 
be destroyed. A sensitive measurement of the spatially and temporally-correlated events 
within the decay series in CsI(Tl) crystals was performed. The measured activities were 
translated to contamination levels of their long-lived parent isotopes in the crystal. If 
secular equilibrium is assumed, the levels of 235 U, 238 U and 232 Th can be derived. 

The identification of the events provides event-by-event background rejection to the 
candidate events for neutrino or Dark Matter- induced interactions. Beyond that, they 
are extremely useful to perform in situ studies of detector performance parameters. 

4.1 Decay Sequences 

The five decay sequences (DS) of interest under studied are [10]: 

DSi) Within the 238 U series, the a-a cascade from 222 Rn — 

222 Rn -> 218 Po + a (Q = 5.59 MeV ; n = 3.82 d) , 
218 Po -> 214 Pb + a (Q = 6.12 MeV ; ri = 3.10 min) , 

the event rate of which gives the concentration of 226 Ra whose half-life (ri) is 1600 y. 
DS 2 ) Within the 238 U series, the (3-a cascade from 214 Bi — 

214 Bi -> 214 Po + z/ e + e" + is (Q = 3.28 MeV ; ri = 19.8 min) , 
214 Po -> 210 Pb + a(Q = 7.83 MeV ; ri = 164 /is) , 

the event rate of which gives the concentration of 226 Ra with ri =1600 y. 

DS 3ai 36) The decay sequence 224 Ra— > 220 Rn^ 216 Po within the 232 Th series can be iden- 
tified by a triple a-a-a signature — 

5.79 MeV ; ri = 3.66 d) , 
6.41 MeV ; ri = 55 s) , 
6.91 MeV ; n = 0.15 s) , 

2 

the event rate of which gives the concentration of 228 Th with ri =1.9 y. 



224 Ra - 


^ 220 


220 Rn - 


^ 216 


2i6 Po _ 


^ 212 



i?n + a (Q = 
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DS 4 ) Within the 232 Th series, there is 64% branching ratio for 212 Bi to decay via a (3-a 
cascade — 

212 Bi -> 212 Po + v e + e~ + 7's (Q = 2.25 MeV ; ri = 60.6 min) , 
2i2p o _^ 208 p6 + a = 8>g5 MeV , . r ^ = Q 3Q ^ ^ 

the event rate of which also gives the concentration of 228 Th. 

DS5 ai 5fc) Though not expected to be dominant, the decay chain of 223 Ra-^ 219 Rn— * 215 Po 
in the 235 U series has a triple a-a-a sequence — 

11.4 d) , 
3.96 s) , 
1.78 ms) , 

the event rate of which gives the concentration of 227 Th with ri=21.8 y. 



™ Ra ^ 219 Rn + 
219p n ^ 215p + 
215p G ^ 211p b + 



a (Q = 5.97 MeV ; n = 

2 

a (Q = 6.94 MeV ; n = 

2 

a(Q = 7.52 MeV ; n = 



4.2 Event Selection 

A total of 40 crystals with a data size of 1725 kg-day taken at the KS Lab was studied. 
Various selection criteria for the decay sequences DS]_ 5 were devised and their efficiencies 
were measured, as summarized in Table 1. With the exception of the fast-cascade in DS4, 
the events were identified as a or e/7 (ID); the delayed events should be within four half- 
lives from the precedent ones (At); the measured energy of the a-events should match that 
of expectations from Quenching Factor measurements (E Q ) while the correlated events 
should have the same Z (AZ=0). Events from DS4 were identified by the convoluted (3-a 
pulse shapes. 

The energy spectrum of the a-events is depicted in Figure 5. The measured energy 
is "quenched" due to recombinations between the ions and electrons such that it is less 
than that of 7/electrons at the same kinetic energy. The quenching factor measurements 
are discussed in details in Section 4.4. The small peak at 6 MeV measured energy is due 
to the energy sum of those (3-a pairs in DS4 whose At are too small to be de-convo luted. 
Apart from this structure, the a-spectrum shows a broad peak at about 2.5 MeV. The 
individual mono-energetic a-emissions within the 238 U and 232 Th series cannot be resolved. 
This poses difficulties and limits the accuracies to radiopurity measurements based only 
on the a-spectra without the timing information. 

The single a rate is 0.0007 kg _1 s _1 or one event per crystal per 11.9 min. This defines 
a time-range under which the identification of correlated a-pairs would be possible, and 
justifies our selection of the five decay sequences for investigations. 
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DS Signatures Decays Selection Background Half-Life (ri) 

Efficiency -to-Signal Nominal Measured 

f a-a 222 Rn-^ 218 Po 

^ 214 Pb 0.93 0.51 3.10 min 3.2±0.2 min 

~2 214 Bi ^214 Po 

^ 210 Pb 0.77 3.2 x 10~ 3 164 fjs 163±8 /is 

!te ~d 224 Ra ^220 Rn 

^ 216 Po 0.86 55 s 54.4±2.4 s 

3b a ^ 212 Pb 0.97 9 x lO" 5 0.15 s 0.141±0.006 s 
~4 ~iTa 2i2 Bi ^212 p Q 

^ 208 Pb 0.37 3 x 10- 5 299 ns 283±37 ns 
"5a" ~a 223 Ra ^219 Rn 

^ 215 Po 0.78 - 3.96 s No Signal 

5b a ^ 2n Pb - - 1.78 ms DAQ Inactive 

Table 1: Summary of the selection efficiencies, background levels and measured half- lives 
of the decay sequences. 



The At distributions for the correlated-events from DSi i2 ,3 a ,3b are shown in Figure 6a- 
d, respectively. Events with delay up to four nominal decay half-lives were sampled, 
except for the DS2 (3-a pairs in Figure 6b, where the maximum delay between triggers is 
500 /j,s constrained by the DAQ settings. In Figure 6a where DSi a-pairs with delay time 
as long as 12 minutes were selected, background due to accidental pairs and those from 
DS,3 a with missing DS^ were statistically subtracted. Selections for the other sequences 
followed from the criteria discussed. The At distributions are compared to exponential 
decay profiles, and the best-fit half-lives are tabulated in Table 1. All measurements are 
in excellent agreement with the standard values. 

The sequence DS 4 has ri of only 300 ns, such that both pulses will appear convoluted 
within the same readout time window of 12.5 /is. An example of such signatures is 
displayed in Figure 7a. The cascade was therefore selected not by timing but by pulse 
shape analysis. The time difference between the f3 and a leading edges is depicted in 
Figure 7b, where the best-fit ri is 283±37 ns, also in good agreement with reference 
values. As shown in Table 1, only 37% of the /3-a events can be resolved. The remaining 
ones are convoluted and contribute to the 6 MeV peak in Figure 5. 

The sequence DSs^s^ is another triple a-cascade. However, the second decay with 
ri=1.78 ms is well inside the DAQ dead time window (starting from 500 /is after the 
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trigger for a duration of about 5 ms), such that this cannot be measured. The first decay, 
with ri=3.96 s, is measurable. The time-difference plot is shown in Figure 8. No evidence 
of an exponential decay was observed, such that only an upper limit to the 235 U series 
can be derived. 



4.3 Radiopurity Levels 

The measured intrinsic radiopurity levels can be presented in several ways and they are 
summarized in Table 2. The measured activity in mBq kg -1 can be converted to radiop- 
urity level in units of "gram per gram (g/g)" of the parent isotopes in CsI(Tl). Assuming 
secular equilibrium, the contamination levels of 238 U and 232 Th can also be derived and 
compared to other measurements. Results were derived from the 1725 kg-day of data 
taken with 40 crystals, and the quoted errors denote the uncertainties of the measured 
mean activities. 



DS Measured Activity Contaminations of Contaminations of 

(mBq kg -1 ) Long-Lived Parents (g/g ) Series''" (g/g) 

238 U : 0.86 ±0.03 x 1(T 12 
238 U : 0.82 ±0.02 x 10~ 12 
232 Th : 2.23 ±0.06 x 10~ 12 
232 Th : 2.3 ±0.1 x 10~ 12 
235 U :< 4.9 x 10~ 14 



1 
2 

3a, 3b 
4 

5a 



0.0107 ±0.0004 
0.0102 ±0.0003 
0.0090 ±0.0002 
0.0061 ± 0.0003 
< 0.003 



226 Ra : 2.92 ±0.11 x 10~ 19 
226 Ra : 2.79 ± 0.07 x 10~ 19 
228 Th : 2.97 ±0.08 x 10~ 22 
228 Th : 3.1 ±0.2 x 10~ 22 
227 Th :< 1.6 x 10- 21 



t assume secular equilibrium 
Table 2: Measured activities and the derived contaminations of the four decay sequences. 



The most accurate measurements were derived from DS 2 and DS 3ai3 fc giving, respec- 
tively, average activities of 0.0102 ± 0.0003 mBq kg" 1 and 0.0090 ± 0.0002 mBq kg" 1 , 
corresponding to contaminations of 2.79 ± 0.07 x 10~ 19 g/g and 2.97 ± 0.08 x 10~ 22 g/g 
for 226 Ra and 228 Th, respectively. Assuming secular equilibrium, the 238 U and 232 Th con- 
centrations in CsI(Tl) are at the 8.2 ± 0.2 x 10~ 13 g/g and 2.23 ± 0.06 x 10~ 12 g/g levels, 
respectively. The scattered-plot for the 238 U and 232 Th levels among the 40 crystals are 
shown in Figure 9. The RMS spreads are, respectively, 80% and 88% for the 226 Ra and 
228 Th levels — or alternatively the 238 U and 232 Th contaminations in the case of secular 
equilibrium. No strong correlations between the two series were observed. 

Comparison with other radiopurity measurements in CsI(Tl) is given in Table 3. The 
measured levels reported in this article are consistent with and improve over independent 
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measurements using high purity germanium detectors on samples of Csl powder. Only 
upper limits of < 7.6 x 1(T 9 g/g and < 7.3 x 1(T 10 g/g could be derived for the 238 U and 
232 Th series, respectively, from the absence of their associated 7-lines. The measurements 
are also consistent with the previous results on CsI(Tl) crystals [8, 9], where the radiopu- 
rity levels were derived by fitting the various a-activities to the measured a-background 
spectra. Reduction of the measurement uncertainties is due to the stringent selection 
criteria for the time-correlated pairs which greatly suppress the background. 
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0.38 


1.0 




[12] 


DS 2 , 3a 


4.0±0.2 


2.4±0.5 




[13] 


DS 2 , 3a 


3.0±0.1 


17±1 


CdW0 4 


[14] 


DS 2 , 3a 


< 0.4 


9.4±0.7 


GSO(Ce) 


[15] 


DS 2 , 3a 


110±20 


2.9±0.3xl0 4 


CeF 3 


[16] 


DS 2 , 3a 


< 5.6 x 10 3 


1.4 x 10 4 


BaF 2 


[17] 


DS 2 , 4 


1.1 x 10 5 


10 5 


CaF 2 (Eu) 


[18] 


DS 2 , 3a 


8920±80 


24.3±0.4 


Liquid 










Scintillator 


[19] 


DS 2>3a 


3.2 x 10" 4 


< 2.5 x 10~ 4 



Table 3: A summary of the measurements of 238 U and 232 Th concentrations in various 
crystal scintillators, assuming secular equilibrium within the series in all cases. Most 
results are due to measurements of the correlated events from the various decay sequences. 
HPGe denotes measurements of 7-peaks with a high-purity germanium detector, while 
4> a represents statistical fits of the purity levels with the measured a-spectra. 

A comparison of the measured intrinsic radiopurity levels of 238 U and 232 Th with those 
in other crystal scintillators as well as in the kilo-ton low-background liquid scintillator 
CTF (Counting Test Facility) detector [19] at Gran Sasso is given in Table 3. The purity of 
the CsI(Tl) crystals with respect to these series is comparable with the best measured for 
Nal(Tl) and CdW0 4 and exceeds those of the other crystals. Nevertheless, contaminations 
in all crystals are much higher than those in liquid scintillator. 
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The sequences (DSi, DS 2 ) and (DS 3a)3b , DS 4 ) originate from the same long-lived parent 
isotope: 226 Ra and 228 Th, respectively. As shown in Table 2 The agreement between their 
measured activities are excellent. Alternatively, the ratios between DS 4 to DS 3a)3 f, give 
the branching ratio of the 212 Bi ^ 212 Po ^ 208 Pb branch within the 232 Th series. The 
measured branching ratio is 68±4%, consistent with the expected value of 64%. 

The absence of signals for DS 5a , as shown in Figure 8, was converted into upper 
limits of 0.003 mBq kg -1 in the measured activity, corresponding to 1.6 x 10 -21 g/g in 
the contaminations of 227 Th. Equilibrium within 235 U series imply an upper limit of 
4.9 x 10~ 14 g/g can be derived for 235 U contaminations in the CsI(Tl) crystals. 

Given in Table 4 is an alternative view to assess the merits of our measurements by 
projected sensitivities in the hypothetical scenario where no correlated events are observed 
and only upper limits can be derived. This was indeed the case for DS 5a . The accidental 
rates are due to uncorrected event-pairs which survive the various selection cuts and fall 
into time windows corresponding to four half-lives. One can derive from the rates the 
limiting sensitivities in mBq kg -1 , as well as the upper limits of the radiopurity levels 
of the decay parents and the two series. The limits represent the levels above which the 
238 U/ 232 Th contaminations of a test sample can be detected by the present methods at 
the measured accidental background rates. The restrictive selection criteria adopted in 
this work (ID+At+AZ+E a ) greatly suppress the accidental background and enhance the 
the detection sensitivities for the correlated decay sequences in DS 2)3a ,3b,4- The achieved 
levels of ~ 10~ 16 g/g are comparable to those of massive liquid scintillator [19] listed in 
Table 3. 

These cascade events were used to provide measurements of several important pa- 
rameters on crystal properties and detector performance. In the following sub-sections, 
we adopted for analysis a subset of the crystals (6 and 13 crystals out of 40 for the 
238 U and 232 Th series, respectively) which have sufficiently large event samples to provide 
statistically meaningful measurements. 



4.4 Quenching Factor Measurements 

The ionization charge densities are much higher for a-particles than those for electrons 
or photons, giving rise to enhanced recombination rates and therefore the light yield is 
reduced. The quenching factor(QF) of a scintillator is the ratio of measurable light yield 
compared to that due to electrons/photons of the same energy: 

QF = ^ (2) 
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DS Accidental Rate Limiting Sensitivity 
(kg-May- 1 ) (mBq kg" 1 ) 


Contaminations of 
Long-Lived 
Parents (g/g ) 


Contaminations of 
Series* (g/g) 


1 0.53 

2 0.0029 


6.2 x lO- 3 

3.3 x 10~ 5 


226 Ra :< 1.7 x 10~ 19 
226 Ra :< 9.0 x 10~ 22 


238 U :< 5.1 x lO" 13 
238 U :< 2.7 x 10~ 15 


Combined for 238 U Series: 






238 U :< 2.7 x 10" 15 


3a,3b 7.5 x 10~ 5 
4 1.6 x 10~ 5 


8.5 x 10~ 7 
1.9 x 10~ 7 


228 Th :< 2.9 x lO" 26 
228 Th :< 9.4 x 10~ 27 


232 Th :< 2.2 x 10~ 16 
232 Th :< 7.0 x lO" 17 


Combined for 232 Th Series: 






232 Th :< 5.3 x lO" 17 


5a 0.26 


0.003 


227 Th :< 1.6 x lO" 21 


235 U :< 4.7 x 10" 14 



assume secular equilibrium 



Table 4: Projected limiting sensitivities and the contamination levels of the five decay 
sequences in the case where no correlated signals are observed. 



where E ee is the measurable "electron-equivalent" energy of the a-particles and Eq is 
Q-value of decays which is also the kinetic energy of the a's. The quenching is energy- 
dependent, and is usually described by two parameters (a,b) in the form [21]: 

Eee = a V E l • (3) 

l + b-E Q 1 ; 

The identified a-events from the correlated pairs provide a background-free sample for 
the measurement of QF. The value of Eq is exactly known and the entire volume of the 
crystal is sampled, as compared to measurements with an external a-source which are 
sensitive only to a localized sites on the crystal surface. 

Tabulated in Table 5 are the measured QF for the five a's from the various cascades 
studied in this work. The a's from DSi, DS2, DS3 ai 3& are isolated events and their light 
yield can be readily measured by integrating the pulses. The energy measurements for 
a's in DS4, however, involve de-convolutions of the (3-a pulses shown in Figure 7a, and 
therefore are subjected to bigger uncertainties. The energy resolutions of the a-events 
(cr(E cc )/E cc ) are also shown, where cr(E ee ) is the RMS of the respective peaks shown in 
Figure 5. 

It can be seen that the quenching behaviour are different for inorganic crystal scin- 
tillators and organic liquid scintillators. The QF is large in crystal scintillators (~ 0.5 
for CsI(Tl) for a-particles in the MeV range) compared to that for liquid scintillator, 
typically at the QF<0.1 range [25]. In addition, the quenched pulses for a-particles have 
faster decay times compared to those from electrons/photons in crystal scintillators, as 
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DS 



Decays E Q (MeV) E ee (MeV) 



QF 



<r(E ee )/E e 



1 

3a 
1 

3a 
3b 
2 
4 



222 Rn 
224 Ra 

218p 

220 Rn 

216p 
214p 
212p 



218p 

220 Rn 

204p b 
216p 

212 Pb 
210p b 

208 Pb 



5.590 
5.789 
6.115 
6.405 
6.907 
7.834 
8.954 



2.76±0.02 
2.92±0.02 
3.09±0.02 
3.33±0.02 
3.67±0.02 
4.33±0.03 
5.08±0.05 



0.494±0.004 
0.504±0.003 
0.505±0.003 
0.520±0.003 
0.532±0.003 
0.553±0.003 
0.568±0.006 



0.076 
0.074 
0.097 
0.069 
0.068 
0.11 
0.14 



Table 5: The measured quenching factors due to a-emissions in the various decay se- 
quences. 



indicated in Figure 4, while the reverse is true for liquid scintillators. The quenching for 
nuclear recoils in CsI(Tl) at the keV energy range is relevant to Dark Matter searches, 
and have been measured (QF~0.2) by various groups [9, 26]. 

The variations of QF with the kinetic energy of the a's are displayed in Figure 10a. 
The measured quenching parameters a=0.25±0.01 and b=0.33±0.02 MeV -1 provide good 
fits to the data. The RMS variation for the QF among the 13 crystals studied is at the 
6.3% level. The variations of ^-5=) with E ee are displayed in Figure 10b. The energy 
resolution is statistically-limited in four of the data points where a(E ec ) oc VEee- The 
worse resolution of the a from 218 Po in DSi is due to its long half-life such that the 
sample is contaminated by accidental events. Energy measurement of the 214 Po-o; in DS 2 
is affected by the sampling procedures of this particular data acquisition system — the 
pedestals for the delayed events within 500 /xs of the primary trigger were not recorded 
and uncertainties were therefore introduced. Measurement of the 212 Po-o; in DS4 involves 
the de-convolution of the (3-a pulse shapes, which gives additional contributions to the 
energy resolution. 



4.5 Spatial Resolution 

The longitudinal "Z-position" of an event can be derived from the difference of the PMT 
signals between the left (Ql) and right (Qr) sides [6, 27]: 

where R is a dimensionless parameter linearly proportional to Z. Data taken with exter- 
nal collimated 7-sources prior to installation provided measurements of the Z-resolutions 



13 



(a z =2 cm) above 400 keV. This approach is limited by the intrinsic spread of the photon 
interaction sites within the crystal due to multiple Compton scatterings. A better in 
situ method is through the study of the measured position difference (AZ) of the cor- 
related (3-ol and a-a pairs which are emitted at the same site. The AZ distributions of 
fl-a from DS 2 and the two a-a pairs from DS 3a 3b for the combined data are depicted in 
Figure lla&b. Both are centered at AZ=0, indicating that the selected pair of events 
were indeed originated at the same site. The resolution a z , as given by the RMS of the 
distributions, are 2.2 cm and 1.3 cm for the f3-a and a-a events, respectively, while the 
RMS spread of a z among the crystal samples are 19% and 24%. 

The a-a samples provide a more accurate description in the studies of the intrinsic 
spatial resolution, since both of the a-pairs are originated at the same site. The /5-events 
in DS 2 , on the other hand, are accompanied by other 7-emissions such that their exact 
vertices are not well-defined. This explains the worse AZ distribution for the f3-a samples. 

4.6 Contamination Gradient with Crystal Growth Axis 

Crystal growth is itself an purification process. As a result, the impurity concentrations 
in the solution increase during the growth process, such that residual contaminations 
inside the crystals are expected to increase with the crystal growth axis. This is verified 
in Figure 12a&b where the Z-position distributions of the event rates of DS 2 and DS 3ai 3b 
are depicted, corresponding respectively to the gradients of the 226 Ra and 228 Th levels. 
The measured gradients of the combined data are 6.8 x 10 -21 g/g-cm for 226 Ra and 

1.6 x 10~ 23 g/g • cm for 228 Th, while the RMS spread among the samples are 33% and 
104%, respectively. This behaviour is in contrast to the uniform Z-position distribution 
of 137 Cs shown in Figure 2b, since the trace contaminations in 137 Cs cannot be separated 
from the natural 133 Cs by any chemical means. 

4.7 Data Acquisition Dead Time 

The event samples can also be used to measure or monitor the data acquisition dead time 
(DT). The fraction DT/(1-DT) is the ratio of recorded double-a events in DS 3a where 
the third one is missing to the triple-a events where the entire sequence DS3 a and DS35 
are reconstructed. The measured value is DT=9±1%, in good agreement to DT=9.5% 
derived by an alternative method using random trigger events [3]. 
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5 Summary and Conclusion 



We report in this article measurements of the intrinsic radiopurity in CsI(Tl). In addi- 
tion to energy measurements and particle identifications, the selection of spatially and 
temporally correlated event-pairs within the 238 U and 232 Th series can greatly suppress 
accidental background and enhance the measurement capabilities. Sensitivity levels com- 
parable to dedicated massive low-background liquid scintillator detector were achieved. In 
particular, the a-pairs correlated by time as long as 10 minutes were successfully identi- 
fied. The CsI(Tl) crystals measured in this work are among the cleanest compared to the 
other crystal scintillators in the 238 U and 232 Th contaminations. The methods and results 
will be valuable references to other low-background experiments where the suppression 
and measurements of the 238 U and 232 Th background are necessary. 

We are grateful to the referee for pointing out a numerical error. This work was sup- 
ported by contracts 92-2112-M-001-057 and 93-2112-M-001-030 from the National Science 
Council, Taiwan, and 19975050 from the National Science Foundation, China. 
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Csl(Tl) Longitudinal View Cross-Sectional View 




Figure 1: Schematics of the CsI(Tl) target configuration where a total of 93 modules 
(186 kg) is installed for the 2003-04 data taking. 
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Figure 2: (a) Energy spectrum of the entire range and (b) reconstructed Z-position 
tribution of the 137 Cs events with 31.3 kg-day of data. 
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Figure 3: Average pulse shapes for a and (5/^ events 
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Figure 4: Event identification with Pulse Shape Discrimination with: (a) the partial 
charge versus total charge method, and (b) the mean time method. The 7 samples are 
from 40 K, while the a events are with electron-equivalence energy of more than 1 MeV. 
A perfect event identification is thus possible. 
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Figure 5: Measured energy spectrum for the a-events. 
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Figure 6: Delay time (At) distribution for (a) a-a events from 222 Rn— > 218 Po— * 214 Pb 
in DSi; (b) (3-a events from 214 Bi^ 214 Po^ 214 Po in DS2; (c) a-a events from 
224 Ra-> 220 Rn^ 216 Po in DS 3a ; (d) a-a events from 220 Rn^ 216 Po^ 212 Pb in DS 36 . 
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Figure 7: (a) A typical (3-a signal of the cascade 212 Bi— >- 212 Pc— ^ 208 Pb (b) The measured 
time difference between the leading edges in the f3 and a pulses. 
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Figure 8: Delay time (At) distribution for candidate a-a events for 223 Ra-^ 219 Rn— > 215 Po 
in DS5. No time-correlated events are identified in such that only an upper limit to the 
decay rates can be set. 
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Figure 9: The correlations of 226 Ra and 228 Th contaminations among the CsI(Tl) crystals. 
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Figure 10: The variations of (a) the measured quenching factors in CsI(Tl) with the kinetic 
energy of the a-particles, and (b) the figure-of-merit cr(E ee )/\/E ee with the measured 
energy E ee . See text for explanations of the worse resolution in three of the data points. 
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Figure 11: The AZ-distributions from a typical CsI(Tl) crystal on (a) (5-a and (b) a-a 
events, from which the Z-resolutions <t z can be derived. 
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Figure 12: The typical Z-position distributions of the radiopurity contaminations from 
(a) 226 Ra within the 238 U series and (b) 228 Th within the 232 Th series. 
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